Quantum interference effects in particle transport through square lattices 
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We study the transport of a quantum particle through square lattices of various sizes by employing 
the tight-binding Hamiltonian from quantum percolation. Input and output semi-infinite chains are 
attached to the lattice either by diagonal point to point contacts or by a busbar connection. We 
find resonant transmission and reflection occuring whenever the incident particle's energy is near 
an eigenvalue of the lattice alone (i.e., the lattice without the chains attached). We also find the 
transmission to be strongly dependent on the way the chains are attached to the lattice. 
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Quantum interference effects are important in the 
transport of particles in mesoscopic systems. Con- 
sider, for example, a particle traversing through a 
square array of quantum dots. Assume the distance 
between dots is close enough so that the particle 
can hop between nearest neighbor dots. Consider- 
ing only the effect of quantum interference, will the 
particle go through the lattice? Classically, the par- 
ticle has a multitude of paths to go from one end of 
the lattice to the other, depending on the size of the 
lattice. Quantum mechanically, however, construc- 
tive or destructive interference can occur because of 
the different path lengths. Thus, the transmission 
of a particle is not assured even when there are clas- 
sically well-defined paths for it to go through the 
lattice. In this work we investigate the effects of 
quantum interference in the transport of a particle 
in discrete and finite square lattices. 

We consider the particle to be governed by 
the tight binding Hamiltonian from quantum 
percolation^^ . This Hamiltonian has the form 



H 






(1) 



where \i) and |j) represent tight binding basis func- 
tions centered on sites i and j, respectively, and 



ytj 



= 1 ii i and j are nearest-neighbors and Vij = 
otherwise. The sum is only over nearest-neighbors. 
In quantum percolation the particle is confined 
to traverse through disordered clusters constructed 
from the methods of percolation theory'^ with some 
occupation probability p. For p < 1 there is dis- 
agreement whether particle states are localized or 
extended. In a review by Mookerjee, et. aL-, they 
concluded that all states are localized and transport 
is dominated by statistically exceptional necklace- 
like resonant states. Daboul, et. al«^, by calcu- 
lating the moments of distances between pairs of 
lattice sites using series expansion methods, found 
evidence of a transition from exponentially local- 
ized to extended or power-law decaying states with 
an energy-dependent occupation probability thresh- 



old p{E). Recent numerical studies of the scaling 
of the conductance g by Haldas, et. alfS,, however, 
found all states to be localized and no indication of 
a localization-delocalization transition. In this work 
we only consider the limiting case p — 1 wherein all 
sites in the lattice are occupied, i.e., all sites are 
available to the particle through nearest-neighbor 
hops, and all particle states are thought to be ex- 
tended. However, we will show that even in this 
limit, the transport through the lattice is very sen- 
sitive to the incident particle's energy, varying from 
complete transmission to complete reflection. 

To determine the transport properties of a parti- 
cle traversing through the square lattice, we attach 
semi-infinite chains to the left and right sides of that 
lattice. Call the left semi-infinite chain the input 
chain and the right semi-infinite chain the output 
chain. The particle is made incident to the lattice 
via the input chain. If the particle goes through the 
lattice, then it exits via the output chain. Following 
the Landauer-Biittiker formalism.'', the conductance 
of the system can then be determined from the re- 
sulting transmission and reflection amplitudes. Be- 
cause of the semi-infinite chains the corresponding 
matrix equation resulting from Eq. (^ is also infi- 
nite. Daboul, et. al»^ recently described a method 
to transform the infinitely-sized Hamiltonian matrix 
in Eq. Q into a reduced matrix, H', that is finite 
and involves only the lattice and its connections to 
the semi- infinite chains using an ansatz. We are im- 
plementing this method in this work. 

There are various ways of attaching the semi- 
infinite chains to the square lattice. In this work we 
consider two ways. One is by point to point contacts 
and the other is by a busbar connection. In point to 
point contacts the input chain is singly attached to 
the top-leftmost site while the output chain is singly 
attached to the bottom-rightmost site of the square 
lattice. In a busbar connection the input chain is 
attached to all the sites in the left side of the lattice 
while the output chain is attached to all the sites in 
the right side of the lattice. 
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FIG. 1: Plot of the transmission coefficient, T, against 
the incident particle's energy, E, for a 10 x 10 lattice 
with point to point contacts to the input and output 
chains. The diamonds (o) are the locations of the doubly 
degenerate eigenvalues of the isolated square lattice. 



Daboul, et. ali^, proposed the following ansatz: 
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where n = 0, 1, 2, . . .. The ■(/;_(„_|.i) represent compo- 
nents of the wavefunction along the input chain and 
the ■!/'+(n+i) represent components along the output 
chain. ^_i and tp^i are for the sites in the input 
and output chain, respectively, that are directly con- 
nected to the lattice. The ansatz restricts solutions 
to Eq. Q in the form of incident and reflected plane 
waves along the input chain and transmitted plane 
waves along the output chain. Because of this ansatz 
the energy of the incident particle is also restricted 
to be within E = —2 and E = 2. The transmission 
and reflection coefRcients can be determined from 
the t and r in Eq. Q by T = |tp and R ^ |rp. 

Once the Hamiltonian matrix in Eq. |^ is reduced 
to H' , the resulting problem can then be cast into 
the form of a linear equation {H' — E)^ = j, where 
7 is solely a function of E. This linear equation can 
then be solved for -0 once E is chosen. We determine 
t and r from tp by numerically solving the above lin- 
ear equation exactly, i.e., from ip = (H' — E)~^j. 
The matrix [H' —E) is sparse and is numerically very 
close to being singular, making the use of standard 
methods such as LU decomposition fail in some in- 
stances. As such, we implement the technique called 
singular value decomposition^ to carefully determine 
the inverse of [H' — E). 

Shown in Fig. ^ is the plot for the transmission 
coefficient against the incident particle's energy for 
a 10 X 10 lattice with point to point contacts to the 
input and output chains. Also shown are the lo- 
cations of the doubly degenerate eigenvalues of the 
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FIG. 2: Sample states of a particle traversing through a 
20 X 20 lattice with point to point contacts to the input 
and output chains, (a) Highly transmitting state with 
E = 0.39. (b) Highly reflecting state with E = 0.41. 



isolated square lattice. An isolated lattice is one 
where the input and output chains are not attached. 
The system is highly transmitting except at some 
values of energy where there are sharp dips and the 
system becomes highly reflecting. Notice that the 
dips occur near the eigenvalues of the isolated lat- 
tice. This phenomenon is analogous to the resonant 
tunneling^ of an incident particle through, for exam- 
ple, a double-barrier system. In this work, however, 
although there is no tunneling involved, we do see 
resonant reflection whenever the energy of the inci- 
dent particle falls near a doubly degenerate eigen- 
value of the isolated lattice. 

Notice as well that there is symmetry between the 
E > Q side and the E < side. The square lattice 
has bipartite symmetry and point to point contacts 
connections preserve that symmetry. Maintaining 
bipartite symmetry can in turn be shown to ensure 
the symmetry in T about E = Q. 

As the size of the isolated lattice is increased the 
number of its associated eigenvalues will also in- 





FIG. 3: Plot of T against E for a 10 x 10 lattice with 
busbar connections to the input and output chains. The 
diamonds (o) are again the locations of the doubly de- 
generate eigenvalues of the isolated square lattice. 



crease. For the lattice with point to point contacts 
to the chains, we also see more dips in the transmis- 
sion coefficient as we increase the size of the lattice. 
These dips are also located near the doubly degener- 
ate eigenvalues of the corresponding isolated square 
lattice. 

Shown in Fig.|2fa) and (b) are sample states that 
are highly transmitting and highly reflecting, respec- 
tively, for a particle traversing through a 20 x 20 
lattice. The lattice is at the xy plane. The input 
chain is attached to the site in the lattice located at 
(1,1). The output chain is attached to the site in 
the lattice located at (20, 20). The z-axis is the ab- 
solute square of the components of the wavefunction 
at each corresponding lattice site, \'\l'{x,y)\ . For 
the highly transmitting state we see a diagonal line 
of non-zero "0 going from the input to the output 
chains. Though this is not always true for all highly 
transmitting states, those with this feature are al- 
ways highly transmitting. In the highly reflecting 
state, on the other hand, we see large fluctuations 
and destructive interference is manifest at the input 
and output sites. 

Shown in Fig. |2| is the transmission T versus the 
incident particle's energy E plot for a 10 x 10 lattice 
with busbar connections to the input and output 
chains. In contrast to the case for point to point 
contacts, the system is mostly reflecting but with 
sharp peaks in transmission at certain values of the 
incident particle's energy. Also shown in Fig. O are 
the locations, as diamonds, of the doubly degener- 
ate eigenvalues of the corresponding isolated 10 x 10 
square lattice. Notice that for the E <Q side the lo- 
cations of the transmission peaks are near the eigen- 
values of the isolated lattice. This is similar to the 
case with point to point contacts but instead of res- 
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FIG. 4: Sample states of a particle traversing through 
a 20 X 20 lattice with busbar connections, (a) Highly 
transmitting state with E = 0.017. (b) Highly reflecting 
state with E = 0.051. 



onant reflection we see resonant transmission. For 
the £' > side, however, some of the peaks do not 
coincide with the locations of the eigenvalues of the 
isolated lattice. 

From the ansatz shown in Eq. |2Jl, the wave vec- 
tor q of the particle is related to its energy by 
E — 2cos(q), where q = 27r/A. For negative en- 
ergies, the particle's wavelength is constrained to be 
within I < A < 4. For positive energies, the wave- 
length should be within either A < | or A > 4. 
Unlike the case for point to point contacts, the lack 
of symmetry between the E < Q and E > Q sides 
of the plot in Fig. |31 indicates the significance of 
the incident particle's wavelength when undergoing 
through a busbar connection. Mathematically, the 
multiple connections of the busbar destroys the bi- 
partite symmetry of the square lattice, and conse- 
quently destroying the symmetry in T about E — Q. 

Let us call those sites at the sides of the lattice 
that are directly connected to the input and output 
chains as belonging to the input and output connec- 



tion boundaries, respectively. Because of the multi- 
ple connections in a busbar, destructive interference 
can occur at the connection boundaries resulting in 
a vanishingly small transmission through the lattice. 
Some of the minima in transmission in Fig. |21 ap- 
pear to be consistent with rules analogous to opti- 
cal interference minima/maxima conditions on the 
boundary. For example, the condition that an inte- 
ger number of wavelengths fit within the boundary of 
a lattice of size L x L, i.e., the condition L — 1 = nX, 
would suggest that certain values of A result in de- 
structive interference. This would include A = 1, 
i.e., E — 2 from E = 2cos(27r/A), for all L, and 
\ = 2 {E = —2) for all odd values of L. In actual- 
ity, completely destructive interference occurs when 
i? = 2 for all L > 2 and also when E = —2 for all 
L ^ 2 and 4. There are also several other minima 
in T that are consistent with this condition. For ex- 
ample, for L = 5, A = 4 (ii^ = 0) also satisfies the 
condition and indeed it is close to a transmission 
minimum. For L = 6, A = 5/n (n = 1,2,3,4), cor- 
responding to E ^ —1.62 and 0.62, also satisfy the 
condition and they are also close to a minima of T. 
In addition, a,t X = 2 {E — —2) we actually observe 
completely constructive interference for L = 2 and 4, 
where L - 1 = (1/2)A and L - 1 = (3/2)A, re- 
spectively, are satisfied. These observations suggest 
strong influences of interference on or near the con- 
nection boundaries on the overall transmission re- 
garding the busbar connection though this boundary 
interference effect is far from providing a satisfactory 
explanation. In fact, since we have a discrete system 
with unit lattice constant rather than a continuous 
slit as in an optical system, it is not clear why A = 1 
actually leads to destructive interference rather than 
the opposite (except for L = 2). Of course, any influ- 
ence of interference along the connection boundary 
must only be a part of the story since interference 



actually occurs throughout the bulk of the system 
(on most of which A is not even well-deflned) and 
since it must also compete with resonant transmis- 
sion and reflection whenever the values of the inci- 
dent particle's energies at the input chain fall near 
the eigenvalues of the isolated cluster. 

Two sample states for a particle traversing 
through a 20 X 20 lattice with busbar connections 
to the input and output chains are shown in Fig. ^ 
The busbars are connected at the y — 1 and y — 20 
sides of the lattice. Shown in Fig. Efa) is a highly 
transmitting state while Fig.^Jb) is a highly reflect- 
ing state. Notice that the difference in amplitudes 
between the states is several orders of magnitudes. 
In Fig. EJb) strong destructive interference occurs in 
such a way that the state ^ nearly vanishes within 
the lattice. 

In conclusion, we find resonant transmission and 
reflection in the transport of a particle through fi- 
nite square lattices whenever the particle's energy is 
near an eigenvalue of the isolated lattice. The way 
the input and output chains are attached to the lat- 
tice influences the transport behavior of the incident 
particle. For point to point contacts the particle 
is mostly transmitting but with transmission dips 
whenever resonance occurs. For busbar connections 
the particle is mostly reflecting with transmission 
peaks whenever resonance also occurs. There are, 
however, peaks in transmission that can not be ac- 
counted for by resonance. These peaks are results 
of interference originating from the multiple connec- 
tions in a busbar. 
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